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The experimental and theoretical evidence relating spin dependence and long range 
confinement is reviewed. One of the simplest confinement ideas, the dynamic 
electric flux tube picture of Buchmiiller, can be exploited to yield a complete 
effective Hamiltonian. In addition to pure Thomas spin-orbit splitting the other 
relativistic corrections are calculated. It is shown that our effective Hamiltonian 
is identical to the mechanical relativistic flux tube model and gives relativistic 
corrections in agreement with QCD predictions. 



1 Introduction 

The linear growth in energy with increasing quark and anti-quark separation 
has long been accepted as a natural prediction of QCDtl On the other hand, 
the quark's orientational energy has been a matter of continuing discussion. 
It is clear that the long range spin interaction cannot be similar to the short 
range gluon exchange since this would imply a large hyperfine splitting in the 
p-wave charmonium states which is at least ail order of magnitude larger than 
observeda About twenty years ago Schnitzeu considered the ratio of p-wave 
X-state masses, 

R = X2 ~ Xl ~ 0.5 (exp.) (1) 
Xi-Xo 

If the interaction potential were of the Lorentz vector type R > 0.8 but it was 
observed that a long range Lorentz scalar interaction would decrease the ratio 
R and successful phenomenological models could be constructed. The improve- 
ment of course is due to the negative spin-orbit energy from scalar confinement 
whereas the vector case is positive. This observation ignited a long tradition of 
assuming scalar confinement. The R argument assumes the spin dependence 
is given by the Breit-Fermi Hamiltonian in the ladder approximation. The 
argument also assumes that virtual coupling to open flavor channels can be 
ignored. It should also be pointed out that the justification of the hypothesis 
of scalar confinement from a fundamental QCD point of view has never been 
very satisfactory. 



'Adapted from talks given at the Quarkonium Physics Workshop, University of Illinois at 
Chicago, June 13—15, 1996, and at the International Conference on Quark Confinement and 
the Hadron Spectrum II, Como, Italy, June 26-29, 1996. 
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The discussion of long range spin dependence was placed on a firm foun- 
dation by the low velocity expansion of the Wilson loop minimal area law 
by Eichten and Feinbergo For one static and one slowly moving quark the 
spin-orbit energy is given in this formalism by 



2m 2 r \ dr dr 

where V is the static potential and V\ is the long range potential. It was then 
observed by Gromefl that there is a further constraint on these potentials, 

dv = dK _ m (3) 

dr dr dr ' 

where Vi is the short range interaction. At large quark separations V% <C V± 
and with linear confinement V — > ar the spin-orbit energy (||J^) becomes 

/r 8 . .~-_i_s.L, (4) 

2m z r 

which is exactly the Thomas spin-orbit energy term. The "Gromes relation" 
(|J) thus provides further evidence for a pure Thomas spin-orbit interaction at 
large separations. Some doubt has been cast on the original derivation in recent 
work by Williams □ We also await a careful study of the B c p-wsve potential 
by lattice simulation. 

A simple and attractive picture relating confinement and spin-dependence 
was proposed by BuchmullerLl In this picture the quarks carry the color field 
around as they move and for large separations the the field collapses to a purely 
electric flux tube. There are two immediate conclusions one can draw in this 
limit. First, since the color magnetic moments of the quarks are decoupled 
there can be no long-range hyperfine or tensor interactions. Second, since 
there is no magnetic field the only spin-orbit interaction is kinematic, known 
as the "Thomas" spin-orbit energy. In this talk I will extend the Buchmuller 
picture considerably to show that this simple concept leads to a compete theory 
and that it is consistent with other results based in more formal applications 
of QCD. 



2 A Review of Buchmiiller's Argument 

Although QCD is the proper theory of quark interaction, it is not unreasonable 
to discuss some questions within the context of electrodynamics. The flux-tube 
field configuration certainly depends crucially on the nonabelian nature of QCD 
while the interactions of the quarks with the fields does not in our case. This 
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is because we assume a flux-tube field in which only one component of the field 
strength (i.e., the radial color electric field) is non-zero in a co-rotating 
frame. We may thus gauge transform F^ v to an abelian subgroup of SU(3). 

We consider a charged particle of spin S and gyromagnetic ratio g. Its 
classical motion in electric and magnetic fields E and B is given by the Thomas 
equationB If the particle is viewed from a frame in which its motion is only 
radial (the co-rotating frame) the fields in this frame (E' and B') are given by 

E = Ti(E'-vixB'), 

B = 7j_(B^ + vj_ x E') . W 
The spin-orbit energj[3 in terms of co-rotating fields is 

ff,o.=--S-f^-B' + ^- Vl xE'V (6) 

where 7J 2 = 1 — v\ and the charge has been absorbed into the fields. For an 
electric flux tube B' = at every point along the flux tube (and also in the 
quark rest frame). The co- rotating electric field is given by 

E' = -ar (7) 

in the flux tube limit. For a slowly moving quark (» « 1) the spin-orbit energy 
(||) becomes 

H„. .~- S-L, (8) 

where L ~ mr x v. We note this is the same as the Gromes result (^). In what 
follows, the validity of the Buchmullcr picture is assumed as is the consequent 
Thomas-type spin-orbit interaction. 

3 Equivalent Laboratory Fields 

For simplicity we will assume here that the anti-quark is infinitely heavy and 
fixed at the origin. A straightforward extension to a general meson would 
consider two collinear flux-tube segments in the CM system with CM at the 
coordinate origin. As in Fig. 1 an arbitrary point on the tube is r, the quark's 
position is rg, and the momentum conjugate to rg is pg = -iVq. 

As we have discussed, the fields in the co-rotating frame in Buchmiiller's 
picture are 

E' = -ar, () 
B' = 0. [ > 
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Figure 1: Flux tube with static Q at the coordinate origin. The moving quark is located 
at tq. In a non-rotating frame the fields are E and B at the point r on the flux tube. 



We now use the Lorentz transformations (||) to find the lab fields E and B 
which (at the same location on the flux tube) are equivalent to the co-rotating 
fields (g). They are 

E = 7_LE' = -a7 ± f, 

B = aj±v±(r x v_|_) . ^ ' 

For a straight flux tube the magnitude of the perpendicular velocity is 

v 

V± = VQ±, (11) 

r Q 

where vqj_ is the component of quark velocity perpendicular to the flux tube. 

We next proceed to define a four vector potential at points along the tube 
in the lab system. The field vector potential is defined by 

B(r,r Q ,p Q ) = V r xA(r,r Qj p Q ) 

r ( r \ . (12) 

= a — v Q ± 7_l — v Q ± (r x vjj , 

where we explicitly denote the spatial dependence of j±(v±) using (|Tl|). We 
choose the gauge A = A±(r)v±, which implies V ■ A = and A parallel to 
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the direction of motion of the tube. The solution in spherical coordinates is 
then 



A(r,r Q ,p Q ) =vj_^i f dr'r' 2 "ft \ —vq± 
r Q r Jo \ r q 

It is easy to verify that B = (f X vj_)^J^(rAj_) is the same as in (p^). 
Similarly, the time component Ao(r, tq,pq) is defined by 



(13) 



E : ■ V,.,i„ = ' ( ^t>Q± j 0. (141 



giving 



(15) 



^o(r,rQ,p Q ) = a ( dr'j±_ ( — v Q ± 
Jo \ r Q 

Both @ and @ can be explicitly integrated. 

4 Interpretation of the Laboratory Four Potential 

We first observe from ( |l3| ) and (|l5|) that for a pure flux tube the "invariant" 
potential is 

U = Aq-v-A, (16) 



U = a 



o 



dr' 7x «) (1 - v'l) = a f dr 1 y/l - vf . (17) 



The resulting Lagrangian L s tring — — U is exactly that of the string actiorE3liil 
and the dynamics of a straight string with a massive spinless quark at one end 
(the other end being fixed) follows from the Lagrangian 



l-v 2 Q -U{v = r Q ,v Q ,p Q ). (18) 

The above Lagrangian reminds us that in a first quantized effective Hamil- 
tonian all fields are to be evaluated at the quark and all derivatives refer to 
quark coordinates. The spinless quark case is straightforward since no deriva- 
tions of fields appear and U is well defined-. JIhe numerical methods of solving 
this problem are by now well establishedE3E3 An example of the global spec- 
troscopy is shown in Fig. 2. In this example the quarks at the ends of the flux 
tube are assumed massless. 

The potentials we have defined in (jl^) and (|l^) are unusual in that they 
depend on the velocity of the quark. This introduces an apparent uncertainty 
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Figure 2: The global Regge spectroscopy r yourtj ing from the numerical solution of ( jl8[ ) with 
massless quarks (m = 0) and two free endsE3o The vertical axis is orbital angular momen- 
tum and the meson mass M is scaled so that the Regge slope is unity. 

between field derivatives and quark derivatives. We therefore define the re- 
quired quark derivative as 

lim V rQ A M (r,r Q ,p Q ) = lim V r A M (r, r Q , p Q ) . (19) 

r— »rQ r— >rq; 

This assumption ensures a smooth transition from field point to quark and 
in particular it guarantees that the Buchmuller electric tube will yield a pure 
Thomas spin-orbit energy when analyzed in the lab frame. The assumption 
( |l9| ) is satisfied if spatial deriviatives are taken at a constant instantaneous 
angular velocity. 

5 Effective Hamiltonian and Reduction 

The one-fermion Bethe-Salpeter equation can, without approximation, be ex- 
pressed as the one-particle Salpeter equation^ 

A+ [a ■ (p — A) + (3m + A - E] A+t/> = , (20) 
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where the usual energy projection operators are defined by 



A± - E ±Hq/2 E , 

E Q = v/p 2 + m 2 , (21) 
Hq = a ■ p + /3m . 

We suppress the subscript on tq and p<g since it is clear that such equations 
must involve quark coordinates. The effective Hamiltonian is then 

H cS = H + A+(A - a ■ A)A+ . (22) 

We have assumed that the potentials have been promoted to Hermitian opera- 
tions by suitable symmetrization. The use of the Salpeter equation avoids pc&j 
sible inconsistencies inherent in the Dirac equation due to the Klein paradoxic 
For the purpose of determining relativistic corrections either equation gives the 
same result. The proof of this statement follows from the expectation identity 

([E±\F(r)])*0 (23) 

valid to leading order in p 2 /m 2 . From this and the standard reduction of the 
Dirac equation we obtain the reduced (semi-relativistic) approximation, 

H eS ~ m + ~ 7t + Ao + —— h H s , . , 

\ 2m 8m J / om z m ^4) 

H S . Q . = JL- ^L.S-VxA).S. 

2m z r ar m 

To explicitly evaluate this reduction we need to expand in powers of v 2 
and evaluate the derivative at the quark coordinate usin g (|l9| ). For low quark 
velocities we find from the general expressions @ and Ijpthat 



av 2 



A (r,r Q , PQ ) ~ ar+-^r\ (25) 

A(r,r Q ,p Q ) ~ ^r 2 v ± . (26) 
3 r Q 

Using consistently the prescription ( |l9| ) to evaluate the field derivatives we 
obtain 

V 2 A ~ 2a/r Q , 

^ * a, (27) 
ar 

VxA ~ o«Qi(f xvi). 
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Finally, using the fact that a heavy quark carries most of the orbital angular 
momentum 



L ~ mrQV Q1 _ , (28) 



(B4h becomes 



i p 2 p 4 \ a aL 2 a , s 

H eS ~[m+?---f- s )+ar+ - - - — — L ■ S . 29 

to oto / 4to r out T Zm z r 



6 The Mechanical Relativistic Flux Tube Model 

The relativistic flux tube fRFT) model is based on the energy-momentum 
Lorentz transformationE3Ej In the co-rotating frame the electric flux tube is 
assumed to have constant energy per unit length equal to the tension a. In a 
frame where the tube rotates with angular velocity uj = vq^/tq the energy 
of an infinitesimal element is then adr/^/\ — uj 2 r and the total energy of the 
tube is 

which is the same as Aq(v) in (|l5|). Similarly, the angular momentum of the 
tube L t is 

L t au [ r dr'r' 2 
Pt = — = — . , s (r x v ± ) , (31) 
r r J VI - u r' 2 

which is the same as A(r) from (p3|). 
In the free Salpeter equationLL-3 

A + [a • p + /3to- M]A + * = 0, (32) 

the confining interaction is introduced by the "covariant tube substitution'0 

= (E, p) —> (E — H t , p — Pt) , (33) 

giving 

A + [a ■ (p - p t ) + 0m + H f t - Et]A+iJj = . (34) 

This is identical to ( |20| ) since H t = Aq and pt = A. We see that the assump- 
tion of pure chromo-electric field in the co-rotating frame is equivalent to the 
mechanical RFT model. 
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7 Conclusions 

Our main result is contained in the effective Hamiltonian ( p(i| ) and its reduced 
form (p9|). The starting point was Buchmiiller's observation that a moving flux 
tube should be pure electric in the co-rotating frame. Although this makes very 
plausible the pure Thomas spin-orbit term and the lack of spin-spin interactions 
at large distance, it did not directly constitute a complete effective Hamilto- 
nian. The present analysis, done in a non-rotating frame, extends Buchmiiller's 
original conclusions to achieve a complete dynamical theory. In addition to the 
usual kinetic energy terms and the static linearly confining energy, there are 
three relativistic corrections which are dependent on the type of interaction. 
These are the last terms in (p^): the Darwin, the L 2 , and the Thomas terms. 
The Darwin term depends on the method of symmetrizing the Hamiltonian 
and thus is somewhat ambiguous. The L 2 and the Thomas spin-orbit_terms 
are exactly what one expects from a more fundamental QCD approachtH 

An additional significant observation is that the original effective Hamil- 
tonian (|2^) before any semi-relativistic approximations were made is exactly 
that of the "mechanical relativistic flux tube model" This model was con- 
structed by considerations of mechanical momentum and energy conservation 
of a flux tube with quarks at its ends. The relativistic treatment of the present 
case will also then have the desired Regge behavior (for massless quarks) of a 
Nambu stringE3 

Finally, we might mention ±hat one can also start with the direct assump- 
tion that B' = on the quarklia This is another way to build in the Thomas 
spin-orbit interaction and by a suitable Lorentz transform method obtain a 
complete dynamical model. In this case the resulting model is relatively easy 
to solve because the Hamiltonian can be explicitly written in terms of canon- 
ical momenta and coordinates. Although it promises to be a useful model it 
does not contain direct reference to the flux tube properties. The relativistic 
corrections turn out to be different than those in ( |2^ ) and the Regge slope is 
the same as from scalar confinement. 
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